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magnetic flux, spin-orbit coupling strength, gate voltage, charging energy, and interaction parame- 
ters (charge and spin velocity and stiffness). 
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I. INTRODUCTION 

Mesoscopic rings represent an important tool for ex- 
perimental and theoretical studies of various phenomena 
which take place on a submicrometer scale. The ring 
geometry allows one to probe many interesting theoreti- 
cal predictions. One of the most exciting phenomena is 
the generation of geometric phases which are manifested 
in the interference patterns of wave packets propagat- 
ing in the ring. Along with the well-known Aharonov- 
Bohm (AB) effect^ which takes place for both spinless 
and spinful particles, the generation of a spin-dependent 
phase is also possible. This effect, sometimes called the 
Aharonov-Casher (AC) effect^, may occur in the trans- 
port of electrons when they are subject to sufficiently 
strong spin-orbit (SO) coupling. The recent fabrica- 
tion of HgTe rings^i made it possible to directly ob- 
serve the AC phase. In earlier experiments with other 
compounds^*^ the signatures of this effect have been 
also detected. 

In order to probe the AC phase it is necessary to have a 
tool for manipulating the strength of the spin-orbit cou- 
pling. This is provided by the gate-voltage dependence^ 
of the Rashba SO coupling^, which serves as a basis for 
a construction of a spin field-effect transistor^. Chang- 
ing the magnetotransport properties of the ring in this 
way, the experimentalists are now able to study the AC 
effecOii. 

Usually the current through a mesoscopic noninter- 
acting ballistic ring is described theoretically by means 
of the Landauer-Biittiker scattering matrix theory^". 
Geometric phases arising due to both magnetic flux 
and SO coupling can be naturally incorporated in this 
formalismiiiSiiiMii^. Effects of electron-electron inter- 
action and charging energy are not taken into account in 
such a consideration. However, they might be important, 
for example, in small quasi-one-dimensional (quasi-lD) 
rings or in arrays of such rings fabricated in very recent 
experiments^*^. 

In the present paper we calculate the linear tunneling 



conductance of the quasi-lD two-terminal ballistic ring 
with Rashba SO coupling threaded by a magnetic fiux. 
The setup is schematically shown in Fig. ^ The spec- 
trum of electrons in the ring is SO-split into two sub- 
bands. We will assume electron densities at which only 
the lowest radial band is occupied. The electron-electron 
interaction inside the ring is modeled by the parameters 
of the Tomonaga-Luttinger liquid (TLL), the leads being 
nonintcracting. Assuming a weak tunneling between the 
leads and the ring, we compute the leading term of the 
Kubo conductance perturbatively expanded in a series 
of tunneling elements. We mostly follow the approach of 
Ref . [3 where a similar problem for spinless fermions was 
considered. We also make use of the bosonization in or- 
der to calculate the required TLL correlation functions. 
However, instead of the Matsubara formalism, we apply 
the Keldysh real-time approach to this quasiequilibrium 
problem (cf. Ref. ITtI) . Such a combination of the Keldysh 
technique and bosonization appears more efficient for a 
derivation of asymptotic results at temperatures lower 
than the mean level spacing of the ring's spectrum. 

After Ref. Tsl it is known that an electron-electron 
interaction strongly renormalizes the height of tunneling 
barriers between the leads and TLL, and therefore at T = 
electron transport is suppressed. At finite temperatures 
T 7^ the linear conductance vanishes as a power law 
of T, while the effective width of a conductance peak 
grows with T ^ 0. In order to ensure the validity of the 
weak-tunneling approximation, in our studies we assume 
a temperature range where the renormalized tunneling 
rates are smaller than the temperature, Fj^j. <C T . On 
the other hand, finite-size effects remain important at 
T <C 1^0, the single-particle level spacing near the Fermi 
level. 

In the temperature regime Ti^r ^ T <C wo the linear 
conductance is represented by a sequence of resonance 
peaks when plotted as a function of gate voltage and/or 
magnetic flux. In our paper we focus on the problem of 
how the distribution of the conductance peaks depends 
on the external parameters (magnetic flux, SO coupling. 
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FIG. 1: The ring threaded by a magnetic flux "1> is weakly 
coupled to the leads through the tunneling barriers ti and tr 
and capacitively coupled to the gate electrode (Vg). 



gate voltage, charging energy) as well as on the parame- 
ters of the Tomonaga-Luttinger interaction. The pertur- 
bative expansion of the linear conductance in tunneling 
elements is known to break down in the resonance posi- 
tions. Finding the poles of the leading term we can es- 
tablish where the conductance peaks are centered. Thus, 
the study of electron transport in the TLL ring provides 
an effective tool of spectroscopy of its many-body states. 
Conceptually this is analogous to the study of the tun- 
neling conductance between the two parallel quantum 
wiresiS which has been realized experimentally^*^. We 
note that a description of a shape of a particular peak is, 
however, a different problem which is usually tackled in 
a somewhat different manner (cf., e.g., Refs. iHli, 
and it will not be addressed in the present context. 

In our paper we extensively discuss the importance of 
the so-called Klein factors and zero modes (topological 
excitations) of the bosonized Hamiltonian24 for the de- 
scription of distribution of conductance peaks. An ac- 
curate account of the Klein factors is necessary due to 
the presence of spin-orbit coupling in the system. The 
zero-mode sector of TLL decouples from its "continuous" 
(bosonic) sector and contains the whole dependence on 
external parametersiSi2^. The latter appear in the topo- 
logical sector after imposing boundary conditions. We 
elaborate on the procedure of averaging the conductance 
over zero modes in the presence of spin-orbit coupling and 
obtain analytically asymptotic results for the peak posi- 
tions at temperatures lower than the mean level spacing. 
We also reexamine the case of spinless fermions reproduc- 
ing the result of Ref. [3 and discuss it in further detail. 

It is worthwhile to note that the relevance of the topo- 
logical modes for a description of mesoscopic phenomena 
in the TLL rings has been already appreciated in various 
contexts, including studies of persistent2& and Josephson 
currents^-- and the study of the AB phase in chiral Lut- 
tinger liquids^^. The structure of the topological sector in 
the presence of SO coupling has been recently discussed 
as well in applications to persistent22i2S and Josephson2£ 



currents. 

The paper is organized as follows. In Sec. ^ we briefly 
outline the construction of the spectrum of the ring with 
SO coupling. In Sec. Ill we summarize the results emerg- 
ing from an application of the Landauer-Biittiker formal- 
ism to this system. They will be further used as a refer- 
ence in the noninteracting limit. In Sec. IV we present a 
derivation of the Kubo formula in the real-time approach. 
Briefly reviewing the bosonization formalism in Sec. V, 
we derive an expression for the dc conductance to be av- 
eraged over zero modes. The procedure of averaging is 
performed in Sec. VI. We discuss the interplay of the 
externally tuned and interaction parameters in the dis- 
tribution of the conductance peaks, especially focusing 
on the modification of the Coulomb blockade due to SO 
coupling. 



II. MESOSCOPIC RINGS WITH RASHBA 
COUPLING: DISPERSION RELATIONS 

The two-dimensional electron gas with Rashba spin- 
orbit coupling is described by the Hamiltonian 

H = ^{pI +pI) + ania^py ~ cjyp^) + V{r), (1) 

where r — ^/jp+y^. The magnetic field B is intro- 
duced in the kinetic momentum p — > p -I- |A via the 
gauge potential A = y, a;,0). The radial potential 
V{r) confining an electron to the ring geometry can be 
modeled, for example, either by singular isotropic har- 
monic oscillator or by concentric hard walls^^. For these 
or any other types of the radial confinement the result- 
ing quasi-one-dimensional spectrum e„cr(fc) is labeled by 
the radial band index n = 0, 1, . . ., by the angular mo- 
mentum hk = . . . ,—h,0,h, . . ., and by the subband index 
(chirality) a — ±. From now on we will put h = 1. 

If the effective ring's width is much smaller than the 
ring's radius, we can neglect the hybridization between 
the radial bands. We also assume electron densities at 
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FIG. 2: The lowest radial band of the quasi-lD mesoscopic 
ring SO-split into two subbands. 
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which only the lowest radial band (n = 0) is occupied. 
Thus, we effectively consider the strictly one-dimensional 
spectrum (see Fig. ^ which has a parabolic shape and is 
SO-split into two subbanda^S: 

27r2 

e<,(fc) = £o<T(fc) = — 77(fc- /c* - crfcfl)2. (2) 

Here L is the ring's perimeter, A:$ = $/$o is a number 
of flux quanta $0 threading the ring, and the parameter 

depends on the Rashba coupling ur. 

Linearizing the spectrum ^ near the Fermi energy, we 
obtain the four branches 

Erjaik) = ujaik - fc° ) = ujo{k - rjkp - fc$ - crfc/j), (4) 



where 6;/^, 7;/.^ = -^l - 2e,/r, ai/r = -^(1 + li/r), 
and fii/r — — li/r) a-re the phenomenological pa- 
rameters describing scattering in a T-shaped (left I or 
right r) junction. The number of flux quanta is given by 
= ^(fc^ — k\), while the quantity kp = + fc") 

corresponds to the Fermi momentum at zero flux. It can 
be replaced hy kp ^ Nq + ^ , where A/z is a difference 
between the chemical potential of the leads and the Fermi 
energy of the ring, and the integer iVo is related to the 
number 2A^o + 1 of electrons in the ring at A/i = 0. Since 
the expression (j^J is periodic in fc^, the integer part oikp 
can be discarded. Thus, the conductance (O actually de- 
pends on the fractional part of ^ . For future references 
we introduce the parameter A:^ = ^ — i. 

In the weak-tunneling limit e//^ ^ 1 the conductance 
© approximately equals 

^ 4e/er sin^ fci?7r cos^ fc$7r 

27r (cos27rfci? — cos27rfc$)2 + -1(6/ + e^)^ sin^ 2-Kkp 

(6) 

As a function fc^ and , it represents a sequence of Breit- 
Wigner resonances. The conductance peaks occur when 
the resonance condition cos27rfei? = cos27rfc$ is fulfilled 
— i.e., at the values of the parameters 

kp + kg, = npi, kp — kij> = n^, (7) 

where npi and are arbitrary integers. We note that 
in the weak-tunneling limit the resonance condition (0 
remains valid for arbitrary angle xi — Xr between the 



specified by 77 = ± (or r] = R, L) and cr = ±. The Fermi 

angular velocity ujq — (^) ^ defines the level spacing 
of the spectrum and fci;- is the Fermi angular momen- 
tum in absence of a magnetic field and SO coupling. 



III. CONDUCTANCE OF THE MESOSCOPIC 
RING: NONINTERACTING ELECTRONS 

Let us consider the conductance of the ring attached 
to the semi-infinite leads (Fig. For noninteracting 
electrons it can be easily found in the framework of the 
scattering matrix theorjtifi. 

It is instructive to consider first the case of spinless 
fermions with the two linearization points k^^^. One 
finds that in the zero-temperature limit and for the angle 
TT between the junctions to the leads the dc conductance 
readsi^^ 



(5) 

I 

junctions, while the shape of Breit-Wigner resonances is 
quite sensitive to the value of xi — Xr- 

It has been demonstrated in Ref. that for electrons 
with nonzero SO coupling and negligible Zeeman split- 
ting the conductance of the mesoscopic ring is given by 
the sum of the two contributions: G{kp,k^ -\- kn) and 
G{kp, k<s> — ks). In other words, the net effect of the SO 
coupling for noninteracting electrons is the generation of 
the different effective fiux values for the different chan- 




FIG. 3: Splitting of the conductance peaks (solid lines) due 
to SO coupling. The dashed lines correspond to ka = 0. 



G{kp, k^) = — 



16e;er sin^ kpir cos^ fc$7r 



271 [—2aiar + (1 + Jijr) cos 2TTkp — 2f3if3r cos 27rfc$]2 -f (1 — "ii^rY sin^ 2T:kp ' 
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nels. Therefore, the pattern of the conductance maxima 
at fcfl ^ is determined by the resonance conditions 

kp + {k^ ±kji) = nR±, hp - [k^±kR) ^ nL±, (8) 

where nfi± and nL± are arbitrary integers. RecaUing 
that effectively fc^ = /cj? — i, we show in Fig. O how the 
arrangement of the conductance peaks is modified by SO 
couphng. 



IV. KUBO FORMULA 

In order to take into account effects of the electron- 
electron interaction on the distribution of conductance 
peaks, we discuss in this section the Kubo formula for 
the linear conductance. Although this expression is very 
standard, we rederive it using the Keldysh formalism. In 
doing this, we pursue two objectives. First, we would 
like to have better control of the approximations used 
(similar to those made in Ref. Il6|) . Second, we would 
like to deduce an expression for the conductance in a real- 
time representation. Its advantage for the ring geometry 
will be discussed in the next section where the calculation 
of time-dependent finite-size TLL correlation functions is 
concerned. 

In the second-quantized formulation the mesoscopic 
ring attached to the leads is described by the Hamilto- 



H = Hi+Hr + H, + Hj 



(9) 



The left/right lead is described by a Fermi- liquid Hamil- 
tonian Hi/r = J dxc]^^{x) - fij ci/r{x), and the 

tunneling term is Ht = J2iA''^i/rcl/^{xi/r)il'{xi/r)+'R.c.]. 



Here Ci/r and ip are the field operators in the leads and 
in the ring, respectively. 

The Hamiltonian of the central part (ring) He , tp] 
can have any interaction term in addition to the kinetic 
term. In our consideration we will model the electron- 
electron interaction in the ring by the Tomonaga- 
Luttinger liquid which includes only forward-scattering 
processes ("density-density" -type interaction). In the 
framework of this model it is also possible to take into 
account the charging effects. They originate from a ca- 
pacitive coupling of the ring to the gate electrode and are 
described by the Hamiltonian Ec{Nring — ^ with 

the charging energy = e^/2Cg. Here Cg is the gate ca- 
pacitance and Nj-ing is the number operator of electrons 
in the ring. 

The linear response of the system to an applied time- 
dependent bias voltage is described by the Kubo formula 
for the ac conductance 
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dt'e-'''(''-'\[ii{t)Jr{t')]), (10) 



where Ii/r{t) = ie[ti/rc]^^{xi/r,t)'4){xi/r,t) - H.c] is a 
current operator in the left/right junction written in the 
Heisenberg representation. 

In the weak-tunneling limit the expression can be 
expanded in a series of Ht- We make use of the real- 
time Keldysh diagrammatic technique, and for t > t' we 
replace ([/, (i), t(i')]) by 



{Ttii{t)ir{t')) - {ftii{t)ir{t')) = 2lllTi{Ttii{t)ir{t')). 



When expressed on the Keldysh contour, it reads 



(11) 



{Ttii{t)Jr{t')) = {fte'S-^"-^'''^''*''li{t)Tte-'^-^"^^'''^'''"^ (12) 



where the operators without carets refer to the interac- (c) (Gf" — Gf)G^, and (d) GfG^. Here G^y^ and 
tion {Ht) representation. Gf, are the momentum-averaged retarded, advanced. 



Expanding CD to the second order in Ht, we obtain ^nd Keldysh functions of the leads in the real-time 



rep- 



dh f dt[{[[Ii{t),HT{tl)],[HT{t[),Irim 
-oo J —oo 

t l-tl 

dh / dt2{[[[Il{t),HT{tl)],HT{t2)lIr{t')]) 
-OO J — oo 

< r dt'2{[m,[HTit'2)AHTit',),Irit')m. 



resentation 

{G^ - G%/r{t) = ^i{{ci,r{t),c\^M}) = ^imSit)"^, 
Gf^.(t)^-([cv.W,4(0)])=-^^-^ 



l/r 

7 

/r 



' °° and vi/^ is the density of states in the left/right lead at 

The next step is to perform averaging over the leads' the Fermi level. 



states. While doing this, we meet the following combina- 



One can straightforwardly prove that the combinations 



tions: (a) (Gf - Gf )(Gf - G^), (b) Gf (G^ - G^), (a) and (b) vanish identically. The combination (c) gives 
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the following contribution to the conductance: 



M 



'int 



dtidt2- 



1(1 - e-*"*2) 



xRe({[[V;i(0), 4m.Mti)]:€{ti + ^2)}), (13) 

where Fj/^ = 2'kvi/ j.\ti / / {Vi / ^.L) and Vj/^ is the volume 
of the left/right lead. 

From Eq. we derive an expression for the dc con- 
ductance (r2 0) at zero temperature 



27r 



M 



dtidt2 



oo J — oo 



X Re{{[[M0),i^jm>Mti)],4iti+t2)})(U) 

Using the operator identities 

{ [C, AlB} + { [C, B],A} = [C, {A, B}] , (15) 
{C, {A, B}} - {A, {C, B}} = [[C, A],B], (16) 

we rewrite Eq. (|14|) . 



G(") = — r,F,L2 / / dtidt2 

J —OQ J —OO 

+Re{[[M0),i^}mAMti),4it2)}])},il7) 
and further express 

{{[[Mo),^}io)lMti)l4it2)}) 
= {{{Mo).4it2)},{Mti),^lm}) 

+ mit2)A^jio),Mm,Mo)]) 



- {{{^jJ{0),{M0lMti)}},m2)}). (18) 

It is obvious that in the noninteracting limit the only 
term {{{'ipi{0),%pl{t2)},{tpr{ti),ipj{0)}}) survives, since 
the other terms vanish due to the fermionic commuta- 
tion relations. We approximate the dc conductance in 
the interacting case by this dominant contribution 

G « — F,F,i2 / / dtidt2 

X Re{{{M0),4it2)},{A{ti),^j{Q)}}). (19) 

Splitting the four-particle correlator, one can recover the 
formula G « {e^ /2TT)TiTrL'^\G^{u; = 0, a;; - Xr)\^ from 
Ref. E where G^( w = 0, Si — Xr) is a zero-frequency 
retarded Green's function for interacting electrons in the 
ring. This approximation physically means that one scat- 
tering event is completed before another takes place. In 
general, the TLL correlation functions of any order can 
be calculated within the bosonization approach, and this 
approximation can be relaxed. 

The combination (d) with GfG^ also gives a finite 
contribution to the conductance, which, however, van- 
ishes in the noninteracting limit as well. Therefore, we 
will neglect it on the same ground as we have just ne- 
glected the subdominant terms in Eqs. (|17|l and H18|l . 



BOSONIZATION 



A. Spinless case 



In order to compute the four-particle correlator H19|l . 
we will make use of the bosonization technique^*. 

Let us consider for simplicity the spinless case. We 
introduce the shorthand notations for the fermionic fields 
V^i = ip{xi) and V'r = i^{xr), where xi and Xr are the 
angle coordinates of the left and right junctions. In the 
following we assume that xi = and Xr = tt. 

In the bosonization the fields /r are represented as a 
sum of the right- (77 = +, or i?) and left- (ry = — , or L) 
moving components. 



and each of 4'i/r,Tj consists of a topological part -Fl/r,?? 



and a bosonic part ip^,. commuting with each other: 



[i^,Vl = 0. 

The bosonic part is given by 



v!;(x) 



1 



/ La 



(21) 



-i'qkx^'\ 



is a small dimensionless cutoff parameter 



where a 

and the operators brjk, 
tion relation [6^^, ^'^'i./l = S., 



bti. satisfy the bosonic commuta- 

'Skk 



The topological part is important for the finite-size 
TLL with periodic boundary conditions. It includes 
Klein factors F^j, zero-mode operators Njj, and the lin- 
earization points /c," (see Fig. |21 assuming fc^ = 0): 



(22) 



The zero-mode operators N,j — take integer values, 
and the following relations are satisfied^: 



[F^,N^'] = (23) 
- 2S,,,, (24) 
{F^,F,y} = {F^fU^Q for f^^rf. (25) 



The bosonized TLL Hamiltonian Htll = He = Hi, + 
Hq consists of a "continuous" (bosonic) Hi, part and a 
topological Hq part which are decoupled from each other. 



Therefore, the factorization oitpi/^^^ into Fi/j.^^^ 
takes place at any time instant: 



and ^f/,^^ 



i'l/rit) - Fi/r,Rm^^,^j,{t) + Fi/,^i(i)V^/,,i(t), (26) 

where the time evolutions of V'i/rT;(^) and Fij^^,^{t) are 
governed by H^ and Hq, respectively. By the same reason 
the statistical averagings in both bosonic and topological 
sectors are independent of each other. 
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The bosonic part of the TLL Hamihonian is given by 

(27) 



a=l,2 k=l 



where v is the so-cahed charge velocity (the renormahza- 
tion of the Fermi velocity vq = -35? )• The operators dak, 
^ofc ('^ = 172) are obtained from 6^^, 6^^, by the canonical 
transformation ljAl|) . The latter depends on the interac- 
tion parameter 1 — ■^{j^ + K), where K is the so-called 
charge stiffness. For repulsive interactions K < 1, while 
in the noninteracting limit K = j = 1 and v = vq. 
The topological part of the TLL Hamiltonian is 



where ao.i = ^{v ± A) and 



AEc V ^ vK 

uja Kvq vo 



(28) 



(29) 



The topological numbers N,) — Nrj — /c,, are shifted by 
kjj — + Skfj_ , where 

_ 4E,{^^CgVg-2No) + 2Afi-u;o 

redefines the linearization points fc° in order to include 
the dependence on A/i and the gate voltage Vg. In the 
basis N ^ Nb. + Nl, J ^ Nr- Nl, the Hamihonian Hq 
acquires the diagonal form 



(31) 



where N ^ N - 2fc^, J ^ J - 2kg,, and fc^ = TYq + (^fc^- 
One can observe that the whole dependence on A/i, Vg, 
and $ is included in the topological sector. 

Using the commutation relations (|23|l we find the time 
evolution of the Klein factors 



where 



Prj = 2aoN^ + 2aiiV_,, = y [£>iV ± A J] . (33) 



The details of the time evolution of the bosonic fields 
are presented in Appendix^ In fact, they are not very 
important for our purpose. We will only exploit the fact 
that the average of the bosonic fields, 



f{t-n) = WR{t)€\m ^ WlW<L(o)), (34) 

is a periodic function of time which can be expanded in 
a Fourier series 



p=0 



■ipujt 



(35) 



with frequency uj = 2^ and real-valued coefficients 
5p(7). Note that the summation in Eq. (|35|l is performed 
only over non-negative integers. 

The real-time periodicity of (?''(t;7) is inherited from 
the spatial periodic boundary conditions. The occurrence 
of the Fourier series H35(l allows us to perform all time 
integrals explicitly. The analysis of the remaining series 
is a much simpler task. 

Let us make yet another approximation in the spirit of 
Ref. 0. In particular, we split the four-particle bosonic 
correlator in (|19|l . neglecting the anomalous averages 
(e.g., {ijj'^il)^)), the left-right mixing (e.g., ('0iV'/{))i and 
the vertex corrections (averages of operators at the same 
spatial point, e.g., {tfiifi'^)) in the bosonic (continuous) 
sector. At the same time, we do not split the topological 
part of the four-particle correlator (unlike has been done 
in Ref. and perform a single averaging of the whole 
over zero modes. 

Implementing this procedure, we obtain 

{M0)4{t2)Mti)i^lm ~ g''*{t2)g\h) (36) 



where (• • • )z.m. implies averaging over zero modes to be 
discussed later. Collecting all contributions, we find 

G K —TiTrL'^ E 9pi{l)9pAl) / / dtidt2 



Pl ,P2— 

(wpi+ao+P,,) 



-iti{upi+aa~Pr,) 



_|_ gi(Nv+^-ri)-^ ^gi?7ti (t»jpi+ao+P,,) _ g-i?)ti(ijpi+ao-P,,)^ 
X ( g^'i-'nt2(i^P2+a.a+P-r,) _ ^irit2(ujp2+aa-P-n) 



Introducing 



9p{l) 



(38) 



and Ari = + A~ , we can cast Eq. H37|l into the form 

^2 



2tt ^ 



(39) 



2ArAl cos{Nr + Nl + 25fc^)7r 



We remark that the alteration of the angle xi—Xr between 
the junctions would only modify the Fourier coefficients 
gp{j) in Eq. (|38|l as well as the relative phase of the 
interference term ~ ArAl in Eq. (|39|) . Meanwhile, the 
poles of A^ in Eq. (|38l) are not sensitive to the value of 
Xl x^. 

In order to treat further the expression (|39|l we need 
to establish an efficient procedure of averaging over zero 
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modes. But first we are going to discuss the modification 
of the conductance caused by the presence of spin 
degrees of freedom and by spin-orbit coupling. 



B. Spinful case 

Performing a similar bosonization procedure in the 
spinful case, we obtain the following expression for the 
dc conductance 



2 

La 



(40) 



The zero-mode operators N^^a with integer eigenvalues 
are shifted to — ~ kj^a- by fc^o- = -I- (5fc^, 
where 

_ 4E,i^^CgVg-mo) + 2Af,-u;o ^^^^ 

The integer Nq = j ^ fc°^ is related to the number 
4A^o + 2 of electrons in the ring when the parameters 



R<7 — fcicr), 



kB,R = j'^a{kiia ±kLa) 



(42) 
(43) 



The coefficients ao,i = ^{i^c ± Ac -f ± A^) and 02,3 = 
^ (i^c i Ac — i^s T As ) are the components of the quadratic 
form of the zero-mode Hamiltonian 



Ho = 



In the basis 



N, 



Jc.s = {Nr+-Nl + )+cj{Nr^-Nl- 
the Hamiltonian (|48|l becomes diagonal: 



VcNi + Ac Jc' + VsNi + XsJi 



and 



C'cNc + v>^cJc] + (J [lysNs+iiXsJs) ,(52) 



(48) 



{Nr+ + Nl+) + <j{Nr^ +Nl-). (49) 



(50) 



(51) 



where Nc,s = N^.s ~ ■ik^^B and Jc^s = J^^ - 4A:$_7i. 

In Eq. H40() the two components ct = ± seem to be 
independent of each other. However, this is not the case, 
and they are, in fact, entangled due to the nontrivial 
procedure of averaging over zero modes. 



VI. AVERAGING OVER ZERO MODES 



equal zero. The parameter ks vanishes in the absence of 
a Zeeman interaction. The parameter 



(44) 



contains the dependence on A/i and Vg. 

Like in the spinless case, it is convenient to introduce 



8gc 

Wo 



, Acs — 



(45) 



and bJr 



and 7c^s = ilir ^ ^c,s), which are 



expressed through the charge and spin velocities Vs^ 
the charge and spin stiffnesses ^ i^s, and the charging 
energy Ec- 

In Eq. (j40|l the rates Vi and Vr remain the same as 
in the spinless case, since we assume that the density of 
states in the leads is spin independent and equals vi for 
each spin component. The spin dependence appears in 
the functions A^a = A'^^ + A~^ , 



A 



g <:jpAhic)gpAhis) ^^^^ 



^ „ t^cPc + WsPs + ao ± Fjjcr ' 



Pr,a = 2aoiV^o. + 2aiiV_^,^ 



(47) 



A. Spinless case 

The typical expression to be averaged over zero modes 
before the time integration has the form [cf. Eq. H37|) ] 



/„ibiN+ib2J\ 

\^ I z.m. 



(53) 



where &i_2 depend linearly on the time arguments ti 2- 
The trace operation is understood as a summation over 
all possible integer values of Nl and Nr. In the basis 
{N, J) we have to sum over either both even {2m, 2n) or 
both odd {2m -I- 1, 2n -I- 1) eigenvalues. Thus, 



-/3ajo fc^i )^ + A(n— fc^ )^] 



m,7i— — 00 



'f3uioli>{m+^-k^f+\(n+^-k^)^] 



m,n— — 00 



6*3 ( 7rfc,,,e ) 6*3 ( 7rfc$,e '^"oa 



'4 ( 7rfc$, e "'^o^ 



(54) 



where 63^4 are the Jacobian theta functions. Some prop- 
erties of these functions are reviewed in Appendix IbI 



The numerator in Eq. (|53|l can be equivalently rewrit- 
ten in the form 



Tr ( e*''i-'^+'''2-'g-/3^fn 



Tr 



where Hq is obtained from Hq by replacing 

ibi 



k,, 



ib2 



(55) 

(56) 
(57) 



At fow temperatures (3~^ <C loq^Ec, the last exponential 
factor in Eq. (|55|) can be discarded. Using Eq. (|B5|) we 
derive the following expression: 



/ ibiN+ib2.J\ 
\^ / z.rn 



+P2(fc^, fc^)e2^''i/('=M + l/2)+2«f<2/(fc.+l/2)^ (58) 



The "sawtooth" function 



V- (-1)" . o 
^ sm znnx 



n=l 



(59) 



has the period 1 and equals f{x) — —x for x G (^515) 
and /(±i) = 0. The functions Pi. 2(^/^1 fc*) a-re deter- 
mined by 



pi(fc^, fc$) 



1 



4 (^irk^ , 



' e3(7rfc^,e-''^/'3'^oS)e3(7rfe.j,e-''^//3"o^) 
_|_ g/3ajo[i>S2(fe^) + Ag2(fc«)] ' ^^'^■^ 

and P2(fc/j,^*) = Pi{kfi -I- ^,fc$ + |). One can observe 
that pi + P2 = 1 • The function 52 (a;) is introduced in 
(lB9l . 

Let us consider the limit of zero temperature, or 
/3 ^ 00. The expression (|58|) becomes exact in this 
limit. Since 61 and &2 are linear in time, we can per- 
form easily all time integrations. Thus, the averaging 
over zero modes effectively results in replacing N 
2/(^^1 + |'^A'i(2)) and J 2/(fc$ -I- i<5$i(2)), where 
Sfii ~ = and (5^2 = i^$2 which assumes further 
summation over the different topological realizations (1 
and 2) of the ground state with the weight factors pi and 
P2- For /3~^ ^ wq the latter approximately equal 

pi(fcp,/c$) = Q{-vg2{ki^) -\gi(k^)), (61) 
P2ikf,,k^) = <d{vg2{k^) + \g2{k^)), (62) 

and play the role of projectors which divide the elemen- 
tary cell (fc$, A:^) e [— ^, 5] X [— i, i] into two areas (topo- 
logical sectors). 

Let us analyze such a partition of the elementary cell 
and consider the (upper right) quadrant defined by < 
fc^ < ^ and < < i. The function g2{x) = x — 1/i 
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FIG. 4: Shift of the conductance peaks (solid lines) due to 
the charging energy. The dashed lines correspond to the non- 
interacting case (Ac = 0). 



for < X < 1/2, and therefore the border between the 
areas of pi and p2 is given by the equation 



iykn + \ki 



I' + X 



(63) 



For repulsive interactions and, moreover, in the presence 
of Ec ^ the relation i>/X> 1 is fulfilled. 

We can establish the borders between the topological 
sectors pi and p2 in the other quadrants by mirroring Eq. 
(|63|l with respect to the k^ axis, fc$ axis, or both. Thus, 
we obtain that the area of the projector pi is the inner 
part of the elementary cell bounded by the hexagon (see 
Fig. EJ. Respectively, the outer part is the area of p2- 

Let us now analyze the conductance in the upper right 
quadrant. In the inner part pi only the zero harmonic 
{p = 0) of becomes divergent near the border line 
(|63|l . In the outer part p2 the zero harmonic {p = 0) of 
A'^ is divergent near the border line H63|l. From both 
sides of the latter the conductance behaves like 



Gcx 



1 



\i>kn + Xkq 



(64) 



This is an expected result as well as the fact that the pole 
of the conductance (condition for the resonant tunneling 
of an electron) matches with a transition pi — > p2 from 
one topological sector to another. 

In order to identify the positions of the conductance 
peaks, it appears sufhcient to consider just the zero har- 
monics p = of the functions A^, because the higher 
ones (p > 1) do not have any poles at all. 

It is instructive to derive the conductance in the non- 
interacting limit u = X = 1. Using the identity 



00 



i-l)P 



sin TTX ^-^ p + X 

p— — 00 



(65) 



9 



one can show that 



The functions H68|l satisfy the identity 



r,r. 



27r 



sin(^+fc*)7r ■ sm{^-k^)n 



Afj, 



16 



(66) 

This resuh can be recovered in the scattering matrix ap- 
proach, if the width of the Breit-Wigner resonance in © 
is neglected. The expression suggests that the reso- 
nant condition can be satisfied at any fc^ by tuning the 
magnetic flux fc$. 

As was discussed in Ref. fjj, the main qualitative fea- 
ture imposed by an electron-electron interaction and/or 
charging energy is the opening of a window at certain 
values of inside which the resonant condition is never 
met. This situation is shown in Fig. 0] and the corre- 
sponding gap value equals Ac ~ ^(1 — ^). 



B. Spinful case 

In order to implement an averaging similar to 153|l in 
the spinful cases, it is necessary to calculate first the par- 
tition function Tr(e~''^''). The trace operation is now un- 
derstood as a summation over all integer values of iV/f+, 
Nl^, Nji-, and N^-. However, in this basis the Hamil- 
tonian (|48() is not diagonal, and we have to use the ba- 
sis H49II . (15011 instead. The summation rules for the latter 
have been formulated, for instance, in Refs. 25 28. 2I Ap- 
plying them, one can find that the partition function is 
proportional to 



16 

E 



X 03 ^^kB^,e 2*^-0 63 TTkiii,e =""0 



, (67) 



where kxi = kx + jSxi {X — n, B, R) and the summa- 
tion is performed over 16 topological sectors. The latter 
are specified by Sxi given in the table 



i 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


16 







2 


2 





2 


2 








1 


3 


3 


1 


3 


3 


1 


1 







2 


2 











2 


2 


1 


3 


3 


1 


1 


1 


3 


3 







2 





2 


2 





2 





1 


3 


1 


3 


3 


1 


3 


1 







2 





2 





2 





2 


1 


3 


1 


3 


1 


3 


1 


3 



One can define 16 functions (i = 1, . . . , 16) 

P^{kx)^Pl{kx^), (68) 
where pi{kx) equals to 

1 



-1 I v^l6 fl3(7rfcn3) 63{-nkis.j) 83(1^^3]) OsjirkRj) 
^ ^ ^J=2 Bsi^k^) SsiTvki.) e^i-Kks) e^iTTkn) 

1 



(69) 



l + g2/3a;oBj(fcx) ' 

and the functions Bj(kx) are introduced in Appendix [51 



(70) 



At low temperatures /3 ^ <C we have an approximate 
relation 



16 

Pi{kx)^X{Q{-B,{kx)). 

J=2 



(71) 



and the functions Pi{kx) become the projectors which 
divide the elementary cell kx & ^] ^ ' ' ' ^ [— 4, in 



the four-dimensional parameter space into 16 topological 
sectors. 

We can now formulate the rule which prescribes how 
to evaluate the average over zero modes in (00)): it is 
necessary to replace 



4/(fcB0, 



Jc 
Js 



4/(fc$,), 
> 'if (km) 



(72) 
(73) 



and to sum over i — 1, . . . , 16 topological realizations of 
the ground state with the weight functions Pi{kx)- 

Once this procedure is implemented, it becomes suffi- 
cient to consider just the zero harmonics (pc ~ Ps — of 
the functions A^^, Eq. (|46|l . for establishing the positions 
of the conductance maxima. In this respect there exists 
a full analogy with the spinless case, and we refer to the 
corresponding discussion in the previous subsection. 

In the framework of the developed formalism it is pos- 
sible to study the influence of the TLL interaction on 
the distribution of the conductance peaks in the pres- 
ence of magnetic flux and SO coupling. The charging 
effects are also naturally incorporated, and the charging 
energy plays a role similar to that of the TLL param- 
eter Vc. They are both combined into Vc [see Eq. (|45|) ]. 
and therefore the effects produced by each of them are 
analogous. Let us then fix = 1 and vary Ec- In Fig.jS] 
we show the elementary cells of the conductance con- 
tour plot in the (fc<i>, k^) plane for k^ — 0.1 and different 
values of the charging energy. The TLL parameters are 
Vc,s = Acs — 1. One can see how the separate effects of 
SO coupling and charging energy (shown in Figs. 13 and 
01 respectively) merge together. 

In experiments the usual tuning parameters are $ and 
Vg. The parameter Vg appears in the theoretical model 
through both kn and fc^. The Rashba coupling constant 
depends on an applied gate voltage^, which is modeled 



by ur = Oi 



,0 



where Vn 



and K > 



is a dimensionless coefficient. The degeneracy point in 
gate voltage at which the Rashba coupling ur vanishes 

is defined by = Vg{aR = 0) = " . Introducing the 
departure from the degeneracy point Awg = Vg — v'^, we 
then express 



l + K'^{AvgY - 1 



(74) 
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FIG. 5: Conductance peaks for fcjj — 0.1 and different values 
of the charging energy 8Ec/loo — 0.0 (upper left), 0.5 (up- 
per right), 1.0 (bottom left), 3.0 (bottom right). The TLL 
parameters axe Vc,s = Xc,s = 1- 



and 



1 + SEc/ujQ 



(75) 



Here w' = -|- ^ -|- TYn - i determines the shift of the 
whole pattern; we may put at will v' = Q. 

In Fig.iniwe demonstrate the influence of the charging 
effects in (fc$,A?;g) plane. The values of Ec and l'c,s^ 
Ac,s are the same as in Fig. O We observe that upon 
enhancement of the charging energy the gap opens near 
Awg = 0. Due to the presence of the gate-dependent SO 
coupling, the pattern of conductance maxima is more 
complicated than that discussed in Ref. Il6l 



VII. DISCUSSION AND CONCLUSION 

In this work we have studied the tunneling conductance 
of a mesoscopic one-dimensional ring attached to two 
Fermi reservoirs. The interaction inside the ring is de- 
scribed by the Tomonaga-Lutinger liquid. The bosoniza- 
tion approach which is usually adopted for the study of 
such model allows us to include the flux and gate- voltage 
dependence as well as the influence of SO coupling. It 
is remarkable that all these externally tuned parameters 
appear in the topological sector of the bosonized theory. 
Therefore, the accurate treatment of zero modes is re- 
quired in order to describe the mesoscopic phenomena at 
low temperatures /?~^ ^ wq- 

Using the Keldysh formalism we have performed the 
calculation of a linear conductance in the limit of weak 
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FIG. 6: Conductance peaks for k = 1 and different values 
of the charging energy 8Ec/ujo ~ 0.0 (upper left), 0.5 (upper 
right), 1.0 (bottom left), 3.0 (bottom right). 



tunneling between the leads and the ring. The real-time 
approach allows one to obtain asymptotic results for the 
distribution of the conductance peaks in the low temper- 
ature limit. Although the perturbative expansion of the 
conductance in the tunneling strength is not very well 
suited for a description of a shape of conductance peaks, 
it is nevertheless quite efficient for establishing their po- 
sitions. The topological origin of the peaks' distribution 
alludes to its robustness upon small modifications of the 
model. 

We have studied the patterns of the conductance max- 
ima at nonzero spin-orbit coupling as a function of mag- 
netic flux and gate voltage. The Tomonaga-Luttinger 
interaction and the charging energy have been seen to 
contribute in analogous way. In rings of reduced size 
the account of charging effects might appear more exper- 
imentally motivated, and therefore we concentrated on 
their study. In particular, we have made a theoretical 
prediction for the distribution of the conductance peaks 
in the presence of both the charging energy and the spin- 
orbit coupling. We observed an interesting interplay be- 
tween both effects. The SO coupling lifts the degeneracy 
of the conduction peaks, and the charging energy opens 
a gap centered at the remaining points of the degeneracy 
in question. The value of this gap is proportional to the 
charging energy. When the latter becomes very large, the 
Rashba effect is less pronounced. The pattern of the con- 
ductance peaks then approaches the form of hexagonal 
honeycombs which is typical to spinless fermions. 

In conclusion, we have described the interplay between 
Coulomb blockade and Aharonov-Bohm and Aharonov- 
Casher effects for the different values of the charging en- 
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ergy, magnetic flux, and spin-orbit coupling, as is mani- In the zero-temperature limit — > cx) the temperature- 
fested in the contour plots of the tunneling conductance. dependent part of 'D{x) can be discarded, and we obtain 
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APPENDIX A: CORRELATION FUNCTIONS OF 
THE TOMONAGA-LUTTINGER MODEL 

Let us consider for simplicity the spinless case. The 
canonical transformation which solves the two-channel 
TLL model is 



dik 

4k 



bRk 



(Al) 



where u± = a/ (7 ± l)/2. 

The explicit form of the time evolution of the bosonized 
fields reads 



1 



'La 



(A2) 



fe=l 



X [uriDikix - ujt) - u^.^D2k{x + u;t)] , (A3) 



where 



ikx 7t 



D2k = e^^^d 



fex if 



—ik: 



2k 



2k- 



(A4) 
(A5) 



Let us consider the correlation function 



X e 



where 
V{x) 

xix) 



It 



1 — cos kx 



k=l 



-dk 



o/Siuk _ ^ 



E 

k=l 



sin kx 



e-"fc = — In 



1 , 1 - e-"-" 
2i'' 1 



(A6) 

,(A7) 
(A8) 



are the periodic functions of x. Obviously, the function 
HA6|I is also periodic in real time, and therefore it can be 
expanded in a Fourier series with the frequency lu. 



V{x) 



■In- 



(l-e-")2 



2ti (1 - e-*^^")(l - e"-")' 
Hence, the correlation function l)A6|l is equal to 



(A9) 



1 - e 



2 ^i{'qx—ijjt-\-icx) j I ]^ ^—i{r}x-\-Lot—ia.) 

For a; = TT we have 



1 / 1-e-" 



La V 1 -|- e 



(AlO) 



(All) 



Its Fourier expansion H35|) is equivalent to the Taylor ex- 
pansion of the analytic function (1 -I- z)~^ for \z\ < 1. 
Therefore the expansion H35|l contains only non-negative 
Fourier harmonics {p > 0) with the coefficients 



r(7)r(p + i) La 
(-i)pr(p + 7) a-y-i 
r(7)r(p-fi) L 



(A12) 



Computation of a correlation function in the spinful 
case is analogous. A new feature arising in this case is 
the double time periodicity of correlation functions with 
frequencies ujc and ujs , which are in general incommensu- 
rate. 



APPENDIX B: JACOBIAN 9 FUNCTIONS 

The Jacobian function 63 is defined by 

00 

6*3(2,9) = 1 + 2^(7"' cos 2n2, (Bl) 

n=l 

and the Jacobian function 64 can be expressed as 

e4{z,q)^93[z+^,qy (B2) 

Both functions are periodic under the shift z — > z -f tt. 
Making a Poisson summation, one can prove that 



^03(712,6--'/'^). (B3) 



In 



^ ^-aik+zf ^ 

fc— — CXD 

Using the expression for the ratio of two ^3 functions 

'3(zi + z2,g) ,^{-ir 9" 



93(2:1 - Z2,q) 



4 — sin 2nz-[ sin 2nz2, 

ri = l ^ 

(B4) 
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one can establish that 



Hm ■ 



'3 7r(a; + 



ib > 



(nx, 



(B5) 



where f{x) is a sawtooth function introduced in (|59|) . as 
weh as that 



hm 



rri2 ' 

4 . 



^2 



, (B6) 



where mi, m2 = 0, 1, 2, 3 and 



xim2 {^) — ^ ^ 



In the Fourier series ljB7|) one can recognize the func- 
tions 



/ \ nib, , 

cos — (4j: + 7712) — COS — (4a:: + TOi) . (B7) 



51(2;) = 310(2;) -.901 (a;) 



- - 

2 1 2 



9 

16' 



g2{x) = 520 (a;) = -.902(2;) 



{x} 



1 

4' 



53(2;) = ff3o(a;) = -503(2;) 

1 1 3 

X > H — < X 

2 14 2 12 



(B8) 
(B9) 
(BIO) 



1 I 1 
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16' 



where {a;} = (a; modi) is a fractional part of x. The 
function {a:} has a period 1 and possesses a property 
{—a;} = 1 — {a;}. One can notice that 33(3;) — gi{—x). 

The other functions gmim2 slso expressed in terms 
of gi, 52, 93- 

512(2^) = -521(2;) =53(2; -1/2), (Bll) 

532(2:) = -523(2;) =51(2; + 1/2), (B12) 

531(2:) = -513(2;) = .92(2: + 1/4). (B13) 
We also define the following functions: 



D 

^2 


= Vcg2{kfj) - 


1 \ ^ ( U \ 1 

r Ac52(.K*j + 


Vsg2{kB) 


4 Asg2[KR), 


^3 


= Vcg2{kfi) - 


r Ac52(.K*j, 






Bi 


= Vsg2\kB) 


+ Xsg2(kR), 






B5 


= i>cg2{kf_,) - 


f Vsg2{kB), 






Be 


= i>cg2{kf_,) - 


f As52(fcfl), 






B7 


= Ac,g2(fc$) 


+ Vsg2{kB), 






Bs 


= Ac52(fc*) 


+ \s92{kR), 






B9 


= Vcgi{kfj) - 


f Ac5i(fc$) 4 


Vsgi{kB) 


4 As5i(fcfl), 


Bio 


= Vcgzik^) - 


f Xcgzik.^,) 4 


Vs93{kB) 


+ Xsgsikn), 


Bu 


= Vcgzik^) - 


f Xcg3ik<s>) 4 


VsgiikB) 


4 As5i(fcfl), 


B12 


= i>cgi{kf_,) - 


f Acgi(fc$) 4 


Vs93{kB) 


+ As.93(fcfl), 


B13 


= Vcgsik^) - 


f Acgi(fc$) 4 


Vs93{kB) 


4 As5i(fcfl), 


Bi4 


= i'cgaik^) - 


f Ac.gi(fc$) 4 


Vsgi{kB) 


+ As.93(fci?), 


B15 


= i>cgi{kf_,) - 


f Ac53(fc*) 4 


Vsg^ikB) 


+ As5i(fcfl), 


BlQ 


= Vcgi{kfj) - 


f Ac53(fc$) 4 


Vsgi(kB) 


+ As53(fcfl)- 
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